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■ Some of the new developments in the theory of heavy ion collisions arc 
\ reviewed. Much of the last progress have been triggered by the high ener- 

[ — ■ gies available at RHIC. In the near future, the LHC will extend the energy 

reach in heavy ions by a factor thirty and give access to new QCD regimes 
0^ , characterized by large densities and temperatures and corresponding mod- 

ified evolution equations. 

O 

■ 1. Introduction 

The advent of collider energies to heavy-ion physics is leading to a pro- 
found change in the field. The new tools available, specially the access to 
the large transverse momentum part of the spectrum, allow for an unprece- 
dented characterization of the high-density state created in such collisions. 
Accordingly, the traditional goal of producing a Quark- Gluon Plasma, the 
r> ' deconfined state of quarks and gluons predicted from QCD, has been en- 

larged to the study of several other mechanisms as thermalization, parton 
distribution functions at very small- x or in-medium evolution of parton 
showers. In general words, a new line is emerging which attempt to study 
how the collective properties of the fundamental interactions appear. For 
that end, large energy densities in extended regions need to be produced 
in contrast with the traditional direction in high-energy physics which at- 
tempt to create the largest possible energy scales in well localized spacial 
regions for new physics to become observable. 

The relevant questions which can be addressed in experiments of heavy 
ion collisions can be (artificially) classified depending on the time scale of 
the relevant phenomena: i) before the collision, the structure of the two 
Lorenz-contracted nuclei is rather different from a typical hadron at the 
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same energy and collective phenomena, nowadays generically known under 
the name of Color Glass Condensate, appear; ii) once the high-density state 
is created, the relevant question is how thermalization and other collective 
mechanisms, as a possible hydrodynamical behavior, appear; hi) finally, 
the properties of the eventually equilibrated medium need to be studied 
by means of some indirect signals. The LHC is in an excellent position to 
address these questions p. 

1.1. Hard Probes provide a general framework 

Reaching collider energies provides the additional tool of the hard part 
of the spectrum, characterized by large virtualities, to be used for these 
studies. A typical hard cross section can be written in the form 

a AB^h = /a(x1) g2) g, Q 2 ) ^ XuX2t Q2) g, D ^ h (z, Q 2 ) , (1) 

with a factorization between the short-distance perturbative cross section 
a(xi, X2, Q 2 ), computable in powers of a s (Q 2 ) and two long-distance terms, 
the proton/nuclear parton distribution functions (PDF), fA(x,Q 2 ), encod- 
ing the partonic structure of the colliding objects, and the fragmentation 
functions (FF), D(z,Q 2 ), describing the hadronization of the parton i into 
a final hadron h. These long-distance terms are non-perturbative objects, 
but whose evolution can be computed in perturbative QCD. These are pre- 
cisely the objects which will be modified in the case that the extension of 
the colliding system interferes with the dynamics, while the short-distance 
part is expected to remain unchanged if the virtuality is large enough. This 
interference between geometry and dynamics makes hard probes perfect 
tools to characterize the medium properties through the modification of the 
long-distance terms in (|T]). 

A conceptually simple example is the case of the J/^> whose production 
cross section is, schematically, 

a hh^J/V = f i ( Xl ,Q2)®f J ( X2 ,Q2)®a i ^ cS \x 1 ,x 2 ,Q 2 )(0([cc] -» J/*)) , 

where now (0([cc] — > J/^)) describes the hadronization of a cc pair in a 
given state (for example a color octet) into a final J /^. This long-distance 
part is expected to be modified in a medium at finite temperature in which 
the deconfinement avoids the bound states to exist [2]. However, this mod- 
ification, being non-perturbative, lacks of good theoretical control and the 
real situation about the .//^-suppression is complicated [3]. 

From the computational point of view, a theoretically simpler case in- 
volves the modification of the evolution of both the parton distribution and 
the fragmentation functions in a dense or finite-temperature medium. This 
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needs of large scales Q 2 in order for the strong coupling to be small and 
perturbative methods to apply. In general, the presence of a dense medium 
is translated into non- linear terms in the evolution equations. 

2. The initial state: the Color Glass Condensate 

The study of the modifications of the nuclear partonic distributions is 
nowadays known under the generic name of Color Glass Condensate. In its 
original formulation [3] it provides a general framework for the whole colli- 
sion, based on an effective theory separating the fast modes in the nuclear 
wave function from the generated slow modes, associated to small- x gluons. 
This small- x gluons have parametrically large (0(l/a s )) densities and can 
be treated as classical fields. The quantum evolution equation of this setup 
is known and, remarkably, can be written in a rather simple form in the 
large- Nq limit [3]- A sophisticated technology has been developed in the 
last decade in this framework which, in addition to describe the structure 
of the incoming wave functions, aims to provide the link to the subsequent 
evolution into a thermal system [0]. Prom a phenomenological point of view, 
the most successful applications of these formalism are the description of 
the multiplicities measured at RHIC [7J the possible presence of the 
predicted geometric scaling in lepton-hadron data [SJ 1101 [S] ; and the sup- 
pression of inclusive particles at forward rapidities at RHIC , which has 
been predicted as a result of small- x quantum evolution ^2] ■ A particularly 
economic description is presented in Fig. ^ Here, the saturation scale is 
obtained from lepton-proton and lepton-nucleus data as Q 2 at oc x~ x A 1 / 3S 
with fitted parameters A = 0.288 and 5 = 0.79. Fig. ^left shows the quality 
of the geometric scaling in lepton-proton and lepton-nucleus |8, data. 
Assuming the same scaling to hold in AA collisions, the multiplicity in the 
central rapidity can be written as |S] 

= N V~s X N^ t . (3) 

where only a total normalization factor Nq is needed once the energy and 
centrality dependences are fixed by lepton-proton and lepton-nucleus data 
respectively. Fig. ^ shows the comparison of this simple formula with 
available data 13 . 

3. The medium modification of jet properties 

A perturbative modification of the fragmentation of the highly virtual 
produced partons, has been proposed as a tool to characterize the medium 
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T h =Q 2 /Q*, 

Fig. 1. Left: Geometric scaling in lepton- proton [H] and lepton- nucleus [B] data. 
Right: Central rapidity multiplicities in pp and AuAu collisions at different cen- 
tralities ^2] an d the corresponding description from Eq. (J3J) 8 . 

produced in heavy ion collisions [Tl]. These effects on high-p^ particles, 
generically known as jet quenching, constitute one of the main experimental 
observations at RHIC so far j^Hj- In the near future, the LHC will extend 
the range in transverse momentum for at least one order of magnitude [Q . In 
the vacuum, a perturbatively produced high-pj parton with virtuality Q 2 ~ 
Pt develops a parton shower by radiating other partons (mainly gluons) 
with decreasing virtuality. This shower stops when a typical hadronic scale 
C(lGeV 2 ) is reached. The resulting jet is an extended object in which 
the collinear and soft emitted gluons are emitted inside a cone around the 
original parton direction. When the high-pt particle is produced in the 
medium created after a heavy ion collision, this parton shower is modified. 
The jets, being extended structures, provide an excellent tool to characterize 
the medium properties at different scales. 

At high enough parton energies, the main mechanism driving the modi- 
fications of high-pt evolution is the medium-induced gluon radiation. As in 
general high-energy processes, the propagation of the partons through the 
medium can be described in terms of Wilson lines averaged in the allowed 
configuration of a medium. Several prescriptions exist for these averages 
and in the multiple soft scattering limit, the saddle point approximation 
of the Wilson lines define a single parameter of the medium, the transport 
coefficient, q, given by the average transverse momentum squared acquired 
by the gluon per mean free path. The typical energy spectrum of radiated 
gluons is softer than in the vacuum although formation time effects pro- 
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vide a typical scale Cj ~ q 1 ^ 3 under which the radiation is suppressed. In 
the same way, the angular distribution is regulated at angles smaller than 
sin 9 ~ (g/a; 3 ) 1 / 4 . As a result of formation time effects (or coherence ef- 
fects) the medium-induced gluon radiation is, hence, finite in the infrared 
or collinear limits. The two main predictions are the suppression of high-pt 
yields due to additional in-medium energy loss and the associated broaden- 
ing of the jet structures. 

3.1. The suppression of high-pt particles 

Although the degree of theoretical refinement of the jet quenching for- 
malisms is still not completely satisfactory, it provides a good description of 
experimental data with the medium-modified fragmentation function com- 
puted as 

D^ h (z, Q 2 ) = P E (e) ® D^ h (z, Q 2 ) (4) 

here, Pe(c) gives the probability of an additional in-medium energy loss and 
is normally computed by assuming a simple Poisson distribution with the 
medium- induced gluon radiation as input |15| I16j . Once the geometry of the 
system is correctly taken into account, a fit to RHIC data (see Fig. |2J) gives 
the value of the time-averaged transport coefficient q ~ 5. ..15 GeV 2 /fm. 
This large value and the corresponding uncertainty is a direct consequence 
of the surface trigger bias effect in inclusive particle suppression measure- 
ments |17l 1181 119j . This is an intrinsic limitation of inclusive measurements 
on the characterization of the medium and on the study of the dynam- 
ics underlying the propagation of highly energetic partons through a dense 
medium. Further constraints can be found by i) measuring different parti- 
cles species, and in particular heavy quarks, as the formalism predicts the 
hierarchy AE g > AE™=° > AEq^°; ii) by directly measuring the induced 
radiation, i.e. by reconstructing the jet structure in a heavy ion collision. 

New data from RHIC on non-photonic electrons attempt to an- 
swer the question on heavy quark in- medium energy loss. These electrons 
are expected to come from the decays of charm and beauty quarks. The 
perturbative description of the relative contribution of both quarks to the 
electron yield is, however, not under good control and the b/c crossing point 
can be as low as 2 GeV or as large as 10 GeV when the usual mass and scale 
uncertainties are taken into account in the approximation. This translates 
into a large uncertainty in the ratio R\a as can be seen in Fig. |3] |22| . The 
description of the experimental data is reasonable within the error bars, 
although not completely satisfactory. A clear distinction between heavy 
mass effects in medium-modified gluon radiation seems only possible with 
a better identification of the c and b contributions and, ideally, by a direct 
measurement of the D and B mesons. 
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Fig. 3. Left: Comparison of the FONLL calculation of single inclusive electrons 
from pp collisions at y / s=200GeV |21j . Right: The nuclear modification factor of 
electrons with the corresponding uncertainty coming from the perturbative bench- 
mark on the relative b/c contribution. Figure from |22|: data from |23j . 



3.2. Jet shapes in opaque media 

The modification of the jet properties is thought to be the most powerful 
tool to investigate the medium properties as well as the dynamics under- 
lying jet-quenching. From an experimental point of view, the callibration 
uncertainty of true jets in the large background environment of a heavy-ion 
collision is the main issue to overcome . The identification of observables 
with small sensitivity to background subtractions is of primary importance 
|24j . From a theoretical point of view, studying medium-modification of jet 
shapes in the most general case is complicated and likely would need a full 
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Monte Carlo implementation of a in- medium parton shower process, which 
is not yet available. Here, we quote a simplified study [2E| which provides a 
quite good phenomenological understanding of the surprising non-gaussian 
shape found in two particle correlations at RHIC j3U] . The main observation 
is that for energies smaller than Co ~ 2(2g) 1 / 3 the typical medium- induced 
radiation angle reaches its maximum value sin# = 1. In this conditions, the 
medium-induced gluon radiation spectrum has the simple form 

d/ med a s C R 1 

■Z™ (5) 



du dkj_ 16-7T 

Taking a typical value q = 10 GeV 2 /fm, Eq. (J5J) is valid for uj < uj 
GeV. Compared to the corresponding spectrum in the vacuum, 

flT v3,c n 1 



dzdk\ 2ir k\ 

the one in © is softer and the typical emission angles larger, producing a 
broadening of the jet signals. 

For practical applications, exclusive distributions, giving the probability 
of one, two... emissions are needed. How to construct such probabilities, 
using Sudakov form factors 



A(t) = exp 



/"* dt' 




7 J to ^ j 





(7) 



is a well known procedure in the vacuum. In j^H] it is proposed to gener- 
alize this procedure by defining eq. (J7J) for the medium by just changing 
dl v&c j dzdk^_ to dl med / dzdk^_ given by Eq. ©• There are convincing motiva- 
tions to make this Ansatz: (i) it provides a clear probabilistic interpretation 
with the right limit to the most used "quenching weights" |151 116j when 
the virtuality is ignored; (ii) the evolution equations in the case of nuclear 
fragmentation functions in DIS can be written as usual DGLAP equations 
with medium-modified splitting functions P(z) — > P(z) + AP(z) |29| . In 
our case Eq. © would correspond to AP(z) with the appropriate factors. 

Two extreme cases were studied in |28| : (1) one of the particles takes 
most of the incoming energy u\ S> <^2 and, hence, 9\ ~ - "J-configuration" ; 
(2) the two particles share equally the available energy, uj\ ~ L02 and 6\ ~ 62 
- "Y-configuration" . 

Case (1). The splitting probability in the laboratory variables $ and 7? 



dP($,z) 



dz d$ 



as ° R EL cos <I>exp l—EL cos 2 $ 1 (8) 



167T 2 ' I IQtt 

r/=0 
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No smearing 
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Fig. 4. Left: The probability of just one splitting © 28 as a function of the lab- 
oratory azimuthal angle A$ for a gluon jet of E- ]ct = 7 GeV. Different medium 
lengths are plotted for the J- and Y-configurations (solid and dotted lines respec- 
tively). Right: Position of the peaks of the A$-distribution and comparison with 
PHENIX data from Ref. [HO]. 




presents a non-trivial angular dependence, which we stress has been found 
by the same perturbative mechanism and parameters describing the inclu- 
sive particle suppression. The distribution found has two maxima whose 
positions are determined by: 



$ max = ±arccosW ^ (9) 
t, La s Cn 



The angular shape in (JHJ) is very similar to the one found experimentally. 
In |28j a simple model is proposed to take into account the different smearing 
effects from the experimental triggering conditions. The results are plotted 
in Fig. 21 for three different medium lengths and Ej et = 7 GeV. To compare 
with the measured centrality dependence on the position of the maxima 30 

I/O _ 

L = N p ^ Tt was taken and plotted in Fig. see |2E] for details. 

The above explanation for the non-gaussian shapes found experimen- 
tally at RHIC has the appeal of presenting a unified description of all high- 
Pt effects in heavy-ion collisions measured at RHIC. This effects could be 
further amplified for jets developing in a flowing medium |2J. However, 
other mechanisms lead to similar away-side correlations. In the most pop- 
ular one, these angular structures correspond to the shock waves produced 
by the highly energetic particle traversing the medium j^T]. In this picture, 
a large amount of the energy lost must be transferred to the medium al- 
most instantaneously and take part in the hydro dynamical evolution. In 
general the released energy excites both sound and dispersive modes, and 
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only the first ones produce the desired cone-signal (in the dispersive mode, 
the energy travels basically collinear with the jet). The energy deposition 
needed for the sound modes to become visible in the spectrum has been 
found to be quite large jHTj. To our knowledge, no attempt has been made 
so far to describe the centrality dependence of the shape of the azimuthal 
correlations in this approach. Given the fact that the two formalisms rely 
on completely different hypotheses, finding experimental observables which 
could distinguish between them is certainly an issue which deserves further 
investigation. 



4. Hydrodynamics and intermediate-pt 

Together with the high-p^, the most relevant data from RHIC so far is 
the measurement of azimuthal anisotropies in the momentum distributions, 
which is considered as the main check of themalization in the medium. 
As originally proposed the hydrodynamical evolution of initial spacial 
anisotropies lead to momentum anisotropies due to the presence of pressure 
gradients. Although not perfectly, the experimental data from RHIC at low 
and moderate-^ can be consistently described by considering the formation 
of a thermalized medium at the very early stages of the collision r < 1 fm, 
which evolves as a perfect fluid |33| - i.e. follows a hydrodynamical evolu- 
tion with negligible viscosity. This observation has lead to the developing 
paradigm of an ideal fluid being created at RHIC [33] • This state would 
be a strongly coupled Quark Gluon Plasma, with properties far from those 
of the asymptotic case of a gas of free quarks and gluons. Whether this 
interpretation is correct or not is a matter of debate at present. 



4-1- Counting the valence quarks of exotic hadrons 

Very interesting effects appear in the intermediate region of 2 < pj_ < 6 
GeV/c. The most spectacular of them is the appearance of valence quark 
number scaling laws for baryons and mesons: (i) Rcp seems to depend 
only on the valence number of the produced particle; (ii) the elliptic flow 
parameter V2 is universal when plotted as V2(p±/n)/n, n being the number 
of valence quarks. The most successful model to describe these features 
is a two component soft+hard model. In this model, the soft spectrum is 
supposed to come from the recombination of quarks in a medium in thermal 
equilibrium jSH]. Thus for a particle with N valence quarks 



dN N 
d 2 P ± dy 



y=0 



C N M ± 



W) 




P± sinh r/± 
f 



k N (P±) ■ (10) 
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A± = 7T/9q is the transverse area of the partonic medium at freeze-out, r 
the hadronization time, M± is the transverse mass of the hadron, Cat is the 
hadron degeneracy factor and 7^ is the fugacity factor for a given quark k. 
The corresponding parameters can be found in [33]. In ^1] the short-hand 
notation for a N-quark hadron 



k N (P±) = K x 



cosh r]± 
f 



N 

£ 

0=1 



m 2 + 



N 2 



has been introduced with Iq and K\ the corresponding Bessel functions. 

The hard part of the spectrum is basically given by Eq. ^ with a 
simplified energy loss formalism, which constitutes an effective way of pa- 
rameterizing the more refined one given in Section |3J In [33] this model 
has been extended to the case of a 4-quark meson to study the sensitiv- 
ity of these observables to make a case for the discovery of exotic states 
in heavy-ion collisions, in particular for the /o(980). In Fig. |SJ Raa for 
different mesons and baryons are plotted and compared with experimental 
data from RHIC. The conclusion from this figure is that if recombination 
is the right mechanism underlying the baryon/meson difference at interme- 
diate p±, heavy-ion collisions are ideal tools to study the quark content of 
different resonances and to find a definitive answer to the structure of these 
exotic states. 




Fig. 5. Starting from below: Raa for ir, Ks, /o(980) as a ss state, p + p, A + A, 
/o(980) as a 4-quark state and S~ + S + . Data from Ref. |37|. 



5. New ideas 

The field of ultrarelativistic heavy-ion collisions has been, arguably, the 
most active field of research and discoveries in QCD in the last years. It 
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is not unexpected that this will continue in the future, in particular within 
the unexplored regions accessible at the LHC. This activity is also reflected 
in the interest shown from other fields, in particular in the relation of ex- 
perimental observables with string theory computable quantities by means 
of the AdS/CFT correspondence. The literature is already very vast, and 
impossible to quote here all the most relevant contributions; let just us 
mention the calculation of the shear viscosity the heavy quark diffu- 
sion coefficient [SHI; the jet quenching transport coefficient jlU]; the presence 
of shock waves the relation with hydrodynamics |42j . 
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